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Abstract 

We consider stimulated pair production employing strong-field QED in a high-intensity laser background. In an infi- 
nite plane wave, we show that light-cone quasi-momentum can only be transferred to the created pair as a multiple of 
the laser frequency, i.e. by a higher harmonic. This translates into discrete resonance conditions providing the support 
of the pair creation probability which becomes a delta-comb. These findings corroborate the usual interpretation of 
multi-photon production of pairs with an effective mass. In a pulse, the momentum transfer is continuous, leading to 
broadening of the resonances and sub-threshold behaviour. The peaks remain visible as long as the number of cycles 
per pulse exceeds unity. The resonance patterns in pulses are analogous to those of a diffraction process based on 
interference of the produced pairs. 



1. Introduction 

Since Sauter's resolution of Klein's paradox it has 
been known that the vacuum is unstable to pair produc- 
tion in the presence of a homogeneous electric field ex- 
ceeding the critical value Eg = m^/e ~ 1.3 x 10^^ 
V/m, when the energy acquired by an electron travers- 
ing a Compton wavelength equals its rest mass |[Tl|2l. 
The view that the critical field was too large to ever be 
produced in a laboratory has recently been challenged 
with the advent of ultra-high power lasers capable of 
reaching i? ^ 10~^Es (3"^. With plasma tools such as 
high-harmonic focussing |5 1 the Sauter-Schwinger limit 
may come within reach during the next decade. 

While the strongest electric fields now available re- 
main below the Schwinger limit, pairs have nevertheless 
been created in the laser experiment SLAC E-144 ||6|. 
Here, the SLAC electron beam was collided with a low 
intensity optical laser to generate high energy ('probe') 
photons. These probe photons then combined with the 
laser to produce electron positron pairs. While nonlin- 
ear intensity effects were unambiguously detected, the 
experiment did not offer a concrete demonstration of ef- 
fects stemming from the electron mass shift fl]. 

The analysis of the SLAC experiment was based on 
the infinite plane wave results of IS] |9] \W\, as the pi- 



Email addresses: theinzl@plymouth .ac.uk (Thomas 
Heinzl), anton . ilderton@physics . umu . se (Anton 
Ilderton), mattias . marklund@physics . umu . se (Mattias 
Marklund) 



cosecond laser pulse contained around 1000 cycles of 
beam. In the next generation of experiments, to be 
performed at facilities such as CLF and ELI ||4|, laser 
pulses will have femtosecond duration, corresponding 
to 0(1) — 0(10) beam cycles. Thus, the use of ultra- 
short pulses now makes it necessary to go beyond the 
approximation of infinite pulse extent. It is well known 
that the pair creation spectra depend very sensitively on 
the details of the laser field, as is apparent from consid- 
ering even the simplest constant field models - for con- 
stant electric fields see |[TT][l2l[l3][l4l[T5l, for magnetic 
fields (HSlinillSlIIl and for crossed fields see lfT8ll20l . 
The impact of pulse shape and duration has recently 
been extensively investigated in the context of vacuum 
pair production lISTl l22l |231 using purely electric fields 
which remain spatially homogenous but become time 
dependent [j24 , 25 1. These investigations have revealed 
a rich substructure in the spectrum of pairs produced in 
electric fields, caused by finite beam duratiorQ 

Here we extend these finite volume studies to 
photon-laser collisions, i.e. "stimulated" pair creation 
lfT2l . using strong-field QED in this particular context 
for the first time. In contrast to the above approaches we 
do not neglect the effects of the magnetic field, which 
are particularly important at high frequencies fy2\, and 
therefore model the laser by plane waves of finite tem- 
poral (longitudinal) extent. Note also that the plane 
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wave obeys Maxwell's equations in vacuum, unlike a 
time dependent electric field. 

We briefly address transverse size effects. One ex- 
pects the following analysis to be valid provided that the 
incoming high-energy photons only probe the near-axis 
region well within the waist size of the high-intensity 
laser With a tight focus (small waist size), as is usually 
required for high intensity, the produced pairs would be 
rather sensitive to the finite transverse extension of the 
laser beam. On the other hand, a weakly focussed laser 
will appear approximately as a plane wave (of finite lon- 
gitudinal extent) to the pairs. In this case, finite trans- 
verse size effects should be minor, and the plane wave 
(or time dependent electric field) model should be a rea- 
sonable approximation. One finds in nonlinear Comp- 
ton scattering (the crossed process of stimulated pair 
production) that in the strongly focussed case where the 
incoming electron beam radius is larger than the laser 
waist, finite size effects generically lead to broad and 
unstructured spectra. For a sufficiently weak focus, on 
the other hand, spectral substructure remains identifi- 
able EH]. 

The paper is organised as follows. We begin with the 
general expressions for pair production in a plane wave. 
We then briefly review the well known results for an in- 
finite plane wave (IPW). Comparing these results with 
pair production in a pulsed plane wave (PPW), we will 
see how the mass shift arises when the momentum trans- 
fer is quantised. We will also show that the behaviour 
of the momentum transfer is intimately related to the 
fine structure of the emission rate, and that pair produc- 
tion in a PPW is essentially a (multiple slit) diffraction 
process. The emission rates in a PPW exhibit below 
threshold behaviour relative to the IPW case. 

2. Plane waves and pair production 

We model the laser by a plane wave of frequency w, 
i.e. a function of the invariant phase k.x (u = |k|), with 
field strength = fj{k.x){k^al — a^^k^),j e {1,2}, 
and a dot is differentiation with respect to k.x. The po- 
larisation vectors obey a*. a-' — —{e^a? /m?)5'''^ , defin- 
ing an invariant, dimensionless amplitude a. The asso- 
ciated gauge potential is A^{k.x) — fj{k.x) a^. 

The total pair creation probability can, in principle, 
be calculated via the (imaginary part of) the polarisa- 
tion tensor in this plane wave background. The com- 
bined space and time dependence makes this a very 
hard problem, and there has been only limited progress 
so far [33!'34l. It seems simpler to directly calculate 
the scattering amplitude for pair production, via strong- 
field QED using Volkov solutions of the Dirac equation 



ll35l l36l . The tree level amplitude for production of an 
e+ pair of momenta and p'^, stimulated by a probe 
photon of polarisation e and momentum k' is 
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where '^p and 5'^^'' denote the outgoing Volkov elec- 
tron and positron, respectively. The Volkov electron is 
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where p is its asymptotic momentum, Up is a free Dirac 
spinor and S is the Hamilton-Jacobi action. 
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Noting that the integrand in ([TJ depends nontrivially 
only on k.x, we introduce light-cone coordinates ~ 
x^ztx^,x^ :— (a;^, x^), and take the laser to propagate 
in the x^ direction so that k.x = ujx^ . The scattering 
amplitude then reduces to a light-cone Fourier integral 
OYCT k.x of M.{k.x) = UpTp^T_p>Vp' e*^^"^^-?'-', with 
the light-cone three momenta p := conserved, 
as follows from (|2]l being an eigenfunction of p. Then 
([TJ becomes 

Sf, = ^S^{p' + p-k')J d(fc.x) e*(>'+y'^'<')'=^-=^7W , 

(4) 

where we have introduced three boost invariant light- 
cone momentum fractions, y, y' and x' for the electron, 
positron and probe photon respectively. They appear 
through the sum 
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which we will see contains all the physics, and which 
we note is Fourier conjugate to the invariant phase k.x. 
We can now look at specific examples. 

3. Pair production in an IPW. 

It is well known that pair production in an IPW is a 
sum over processes supported on 



(6) 



Here n is the number of laser photons taken from the 
beam and is the average, or 'quasi', momentum of 
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the electron over a laser period which results classically 
from its rapid quiver motion in the laser field [37 1- The 
quasi-momentum may be calculated from the classical 
orbit, and in an IPW is 



9m =Pp 
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(7) 



with flQ the dimensionless intensity parameter, often 
written qq = eE/com in the lab frame, and which 
equals a (a/2) for circularly (linearly) polarised waves 
||38]| . Squaring the quasi-momentum gives the cele- 
brated electron mass shift, — m'^{l + a§) —: 



||39]| . Hence, (|6]), summed over n, describes the multi- 
photon production of a heavy pair with rest masses . 
It follows from squaring ^ that a minimum number of 
laser photons is required, namely the first integer bigger 
than Sg = 2m1/k.k' . Hence, the effective mass blue- 
shifts the energy thresholcj^from 2m^/fc.fc'. 

To explain how these structures (including the mean- 
ing of a number of photons coming from a classical 
background field) arise, we must discuss the relevance 
of the quasi-momentum for the scattering amplitude Q. 
Recall that S — p.x + Ip and Ip depends on k.x = lux". 
For an oscillatory field the integrand of Ip decomposes 
into a constant average (or "light-cone zero mode" |41 J) 
and a periodic, oscillating fluctuation Sip. The average 
over a cycle is the longitudinal component of the quasi- 
momentum, so we write 



S — p.x + {q_ — p_)x + Sip 



(8) 



If we use a bar to denote quasi-momentum fractions, i.e. 
y = q_/k^ instead of y = p_/k_, then Q becomes 

(9) 

where M — Ai but with / — ^ SI. Now, we have the 
Fourier integral of a purely oscillatory, periodic function 
M, which may of course be represented as a Fourier 
series, M(k.x) = ^„ Mn i 



-ink.x 



Plugging this into 



(|9]l we finally obtain 

Sf, = -^,53 (p' + p - k') 5] '5(y + y' - x' - n) . 

~ n 

(10) 

This rather neatly expresses the important and general 
result that the scattering amplitude is a Ad delta comb, 
irrespective of the detailed functional dependence of 



-In nonlinear Compton scattering, tlie mass shift leads to a red- 
shift of the Compton edge 1401 . 



the IPW on k.x (as long as it is periodic). The argu- 
ment of the comb, y + y' — x', contains the longitudinal 
momentum transfer from the probe photon to the pair, 
(jSjl, and also the average transfer from the laser to the 
pair (jsjl which shifts p_ to q^ (p'_ to q'_). Hence, the 
delta comb quantises the total longitudinal momentum 
transfer to the pair, forcing it to take integer values via 
y + y' — x' = n E Z, which is equivalent to 



q_ + q'_ — k'_ ~ nk_ . 



(11) 



Since the IPW quasi-momenta (j7|l differ from the 
asymptotic momenta only in the minus component (and 
since p is conserved), it follows that ( fTO] ) has support on 
precisely the sum of delta functions enforcing (|6]l. So, 
the quantisation of longitudinal momentum transfer has 
the appearance of an integer number of photons being 
absorbed from the laser and used to create heavy par- 
ticles of rest mass . In an IPW, therefore, the mass 
shift and photon number are rather closely related. Note 
that our interpretation, and equation (111, are consistent 



with the quantum optics law that photon number and 
phase are conjugate variables. 

The emission rate, calculated from inherits the 

10 1 with Mn |A/„P, and so 
The support of the delta 



delta comb structure of 



becomes an incoherent sum. 
function in ( fTO] ) or, equivalently, light-cone momentum 
conservation ( [TT| , may hence be viewed as (idealised) 
resonance conditions leading to peaks of zero width, for 
instance, in the triple differential rate, i.e. the rate as a 
function of p'^ and n. Such resonances have been pre- 
dicted before in discussions of purely time-dependent 
electric fields Ell |43l HH |45l |46| which may be viewed 
as simple models for counter-propagating laser fields. A 
more realistic version of the latter (taking into account 
the magnetic field) has recently been studied by numeri- 
cally solving the Dirac equation. A resonance condition 
very much akin to ( [TT] ) has been found |32|. Clearly, 
one expects modifications such as line broadening in a 
more realistic situation, in particular upon taking into 
account finite pulse duration which is the subject of the 
following sections. 

4. Pair production in a pulse. 

A pulsed plane wave is described by a field strength 
F^y — F^j_^{k.x) that goes to zero for sufficiently large 
modulus \k.x\. Typically, it will contain a finite number 
of cycles, N, modulated by a smoothly vanishing enve- 
lope function. For the purposes of this paper we take our 
fields to vanish outside of k.x e P = [0, 27r7V]. One ex- 
pects that in the limit of long pulses, i.e. large iV ^ 1, 
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the physics in a PPW approaches that of the IPW. This 
was the conclusion of Kibble (4T\ who, in the sixties, 
ended a controversy on the unphysical nature of IPWs 
by arguing that they were sufficient to describe the long 
pulses which were then state-of-the-art: the width Acj 
of these pulses in frequency spac^ obeyed Aui ^ ui, 
which is equivalent to ^ 1. However, modern high 
intensity lasers typically have pulses that contain only 
few cycles, N = 0{1) such that an IPW can only yield a 
rather crude idea (at best) of the ongoing physics in this 
case. In particular, the appearance and interpretation of 
quasi-momenta and effective mass as averaging effects 
seem to become questionable. The extent to which these 
ideas can be applied to pair production in pulses will be 
examined in this and the next section. 

The simplest "pulse" is obtained by retaining period- 
icity within P, i.e. truncating the IPW fields. (Edge ef- 
fects may be present if the fields do not vanish smoothly, 
though these should become negligible in the large- 
volume limit, i.e. with increasing N.) Thus, we assume 
for now that F^i^ is regularly oscillating for < a; < 
27riV, corresponding to a wave train of N identical cy- 
cles |'48l. We will drop this slightly over-simplistic as- 
sumption and add smooth envelopes in due course. 

For a PPW the representation Q of the S-matrix el- 
ement remains valid, as the only prerequisite there was 
dependence of the background on k.x. Quasi momen- 
tum fractions then appear through the decomposition 
(|8]l. For the cutoff pulses in question the quasi-momenta 
do not necessarily match those of the IPW, but the mass 
shift typically does. Accordingly, we end up with (j9]), 
but with M no longer being strictly periodic. There- 
fore, performing the phase integral no longer yields a 
delta comb as in (|TOi. Reducing the integral to a sum 



over single cycles, one may show that Q takes the form 



S 



1 , , sin Nt:z ~ 

= — p' + p-k' ^ Mo, (12) 
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where Mq is the contribution from a single cycle, giv- 
ing an envelope, and we abbreviate the longitudinal mo- 
mentum transfer as 

z = y + y'-x', (13) 

from here on. Fourier expanding the single cycle con- 
tribution to aid the interpretation we find 



Mo = 2^ Mq sine 7r(z - 



(14) 



Squaring Sfi, there are no cancellations between terms 
as in the IPW, and so the emission rate is a coherent 
sum, containing interference terms. We recognise the 



product of the sine function in ( 14i and the rapidly os 



cillating ratio of sines in ( 12 1 as the root intensity distri 



bution of light scattering through a finite aperture con- 
taining N slits. Thus, pair production is essentially a 
diffraction process: finite beam duration corresponds to 
a finite aperture, while the N cycles of the pulse corre- 
spond to N slits. The emission spectrum will therefore 
take the form of a diffraction pattern, deriving from the 
interference not of light but of pairs produced at dif- 
ferent points in the beam B9ll . (Compare the proposal 
in 1 50 1 to perform the double slit experiment by polar- 
ising patches of vacuum to create a grating.) Interfer- 



ence is seen through the ratio of sines ( 12i, which has 



maxima when the momentum fraction transfer is an in- 
teger, z G Z, and N — 2 additional surrounding maxima 
of smaller amplitude, resulting in the substructure dis- 
cussed in the introduction (see also the figures below). 

As expected, the peaks signal the persistence of the 
IPW resonances with the previous delta comb now 
acquiring finite widths. The longitudinal momentum 
(fraction) transfer is now continuous, but peaks are still 
located at z = n e Z. The resonance condition 
thus remains valid and one may loosely think of z as 
something like a continuous photon number. As z in- 



^Recall that the Fourier transform of a monochromatic IPW is a 
delta function, corresponding to zero width. 



volves quasi-momentum fractions, see ( 13 1, the peaks 
still look like the production of mass m*, rather than m, 
pairs. However, the continuous nature of z together with 
the line broadening phenomenon implies that the sharp 
IPW cutoff 2m1/k.k' becomes washed out as well: in 
other words, there is sub-threshold behaviour! This 
is completely consistent with Kibble's reasoning since 
a PPW contains higher frequency components than a 
monochromatic IPW which obviously lower the thresh- 
old. In fact, the reduction is quite significant: the thresh- 
old is reduced to 2rr? jk.k' , independent of the geome- 
try of the pulse. 

The threshold reduction can be observed in the de- 
tection of pairs at lower energies than the minimum im- 
plied by (j6]l, i.e. with energy lower than is required to 
produce the effective mass or, equivalently, reach the 
first resonant momentum transfer To see this, and the 
preceding, effects explicitly we turn now to the discus- 
sion of examples 

5. Examples 

In what follows, we will study PPWs representing a 
slight generalisation of a model recently suggested by 
Mackenroth et al. |29|. We consider a linear polarised 
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wave by setting /2 ~ 0, and consider the family of 
gauge potentials 



. - , sin^ f TTi^T I smik.x) k.x G P 
fi{k.x) = { \2N J ^ ' 



10 







otherwise 



(15) 



where the envelope is characterised by an integer K. 



5.1. Example 1: Regular wave train (K — 0) 

We begin with K = corresponding to a finite wave 
train of N cycles |48| i.e. a "mutilated" IPW sharply 
cut off outside of P. The potential ( [T5] l is then con- 
tinuous, while the field strength is fi{k.x) = cos k.x 
inside P, zero outside, and therefore has hard (step 
function) edges. These are unphysical, but we will 
see how their effects come to disappear in the large- 
limit. For this example, the quasi-momenta are as in 
(|7]i, and the mass shift remains the same as in an IPW: 

= (1 + a^/2)rrp- = (1 + al)m^ . 

Fig.[T]gives the triple differential rate for pair produc- 
tion in = 1,2 and 4 cycles of the beam, plotted as a 
function of \p'^\/m for fixed transverse angle = 
and half maximum light-cone component p'^ — aj'/2. 
In accord with the discussion of focussing in the intro- 
duction, the intensity is not too high (this is also useful 
for quick numerical calculation of the rates). Vertical 
(black) lines are the four lowest terms of the IPW delta 
comb (the remaining contributions stand to the right). 
There is a rich structure even in the single cycle contri- 
bution which must be dominated by finite pulse dura- 
tion, and in this case edge, effects. As we go beyond 
one cycle (for which the ratio of sines in ( 12 1 is unity) 



the shape of the differential rate changes significantly. 
When A^ > 1, we have interference effects which are 
analogous to those appearing when one goes from sin- 
gle to multiple slit diffraction: resonance peaks, centred 
on the conservation of quasi-momentum as in (|6]l, with 
N — 2 surrounding subpeaks. As is clear from the prop- 
erties of the sin / sin factor in the rates (and as shown 
in the plots), the strong peaks dominate in the limit of 
a large number of cycles, exactly reproducing the IPW 
delta comb. 

We also observe the threshold reduction, through the 
production of pairs with energy lower than that allowed 
in the IPW: this is seen as the non-zero signal to the 
left of the first possible IPW resonance at \p'^\/m = 1. 
Thus, this simple example contains all the promised fea- 
tures. Let us now move on to more physical examples 
without edge effects. 



Ip'-lI/hi 




Figure 1: Triple differential rate, units of e'^m'^n-./ /32n'^uj' , as a 
function of transverse e+ momentum. Linear polarisation, ao = 2, 
cj' = 250 GeV. A'^ = 1, 2 and 4 cycles of the laser (descending), 
corresponding to 4 fs, 8 fs and 16 fs pulses respectively. Black (ver- 
tical) lines are the IPW delta comb (quasi-momentum conservation 
resonances). The non-zero signal to the left of Ip'^l/m = 1 is pro- 
duction below threshold. 



5.2. Example 2: Pulses with smooth envelope (K > 1) 

For our second example we add a smooth envelope 
function. It should be expected from results on vacuum 
pair production 11211 [51 ] that the resulting spectra will 
be highly sensitive to the details of this envelope, i.e. 
for our example family 



15 1, on the value of K. We 



will see, however, that all the features associated with 
the diffraction pattern and predicted above continue to 
appear in the rates. Following [29l we adopt K = A 
hence adding a sin"' envelope to the previous wave train 
{K — 0). The field strength F^^ is now smoothly van- 
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ishing at the edges of the pulse. The rate for this pulse is 
shown in Fig. |2] Even though the pulse no longer con- 
sists of identical cycles, all the features of the diffraction 
rate remain, (i) There is sub-threshold behaviour rela- 
tive to the IPW results, (ii) The form of the rate changes 
dramatically above one cycle of the beam, with rapidly 
oscillating substructure setting in. (iii) This leads to 
a series of strong peaks. The envelope in ([TSj natu- 
rally smooths out many features of the spectrum, and 
the fine structure of the diffraction pattern sets in a little 
slower than in our previous example, so we have plotted 
1, 4 and 8 cycles of the laser. Reassuringly, our inter- 
pretation of pair production in pulsed plane waves as a 
diffraction-like process, and the resulting features of the 
positron emission spectra, are shown to hold for realis- 
tic pulses as well. For additional confirmation of this, 
we display in Fig. [3] the eight-cycle rate for the pulse 
with the K ~ 2 envelope. The field strength is again 
smooth. The sub-threshold behaviour and interference 
pattern remain, but clearly differ from those in Fig. [2] 
This illustrates both the expected sensitivity to the pulse 
envelope, and the robustness of our results. 

We remark that in arbitrary pulses, there is no single 
natural definition of an average from which to construct 
the quasi-momentum and mass-shift. Assuming iden- 
tical cycles as we did before, the quasi-momenta are 
given by a cycle average, but do not in general agree 
with their IPW counterparts (though the mass-shift typ- 
ically does). Indeed, we notice in the above plots that 
the strong peaks of the rates do not exactly track the 
IPW resonances as in the K ^ example, which is 
an indication of (physical) finite size effects affecting 
the quasi-momentum. However, our diffraction patterns 
offers us a potential way to reconstruct the mass shift. 
This is possible because any separation into an 'aver- 
age' and 'fluctuation' of the field will lead to a separa- 
tion of the Volkov exponent along the lines of (j8]l. This 
results in an expression for z as a function of |p^|/m 
(at each angle transverse angle (p) which is just the reso- 
nance condition for the given pulse. Hence, by locating 
the positions of the resonance peaks in the production 
rates, one could in principle reconstruct this function, 
and from it the quasi-momenta and the mass shift. We 
hope to comment further on this idea, and the mass shift 
in general, in a future paper. 

6. Conclusions. 

We have given a new interpretation and understand- 
ing of the electron mass shift in strong-field pair pro- 
duction. We have seen that pair production in a plane 
wave of finite temporal extent is a diffractive process. 




Ip'J/m 



Ip'J/m 
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Figure 2: Triple differential cross section for the pulse \15) , K = 4, 
field amplitude a = 2\/2 to match the intensity parameter above. 
N = 1, A and 8 cycles of the laser (descending), corresponding to 
4 fs, 16 fs and 32 fs pulses respectively. The diffraction pattern re- 
mains. Sub-threshold behaviour remains. 




Figure 3: Triple differential cross section for the pulse \15\ with K = 
2. Other parameters as in Fig. |2] and 32 fs duration. The rate again 
exhibits both the diffraction pattern and sub-threshold behaviour. 
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The positron emission spectrum can be interpreted as 
an interference pattern, with a rich substructure. The 
rate exhibits resonant behaviour when the laser aver- 
aged Hght-cone momentum transfer to the pair is a mul- 
tiple of the laser frequency. This resonance condition 
is equivalent to the momentum conservation of a multi 
laser-photon (higher harmonic) process which creates a 
pair of effective mass m*. However, the rates are not 
completely suppressed away from these conditions, as 
in the IPW case, and thus the pulsed rates exhibit line 
broadening and significant below-threshold behaviour, 
with the electron rest mass, rather than the effective 
mass, setting the scale. For a large number of cycles 
the diffraction pattern resembles the delta comb of the 
IPW limit, centred on the resonant values of momentum 
transfer. This corresponds to the loss of the coherent 
(interference) terms in the pair production rates. 

Phenomenologically, pair creation certainly requires 
high energy probes to stimulate the process. With re- 
gards to the laser parameters, it seems reasonable to 
work at moderate ao = 0(1) which, in turn, allows for 
longer pulse duration and therefore a large number of 
cycles per pulse, TV ^ 1. In this regime, finite pulse du- 
ration effects should be under control. This will require 
some fine tuning for which the present results should 
provide a solid basis. 

Above, we considered the triple differential rates, 
which, unless the pair yield is very high, would be chal- 
lenging to measure experimentally. Our focus was on 
these rates because they are the simplest to plot and un- 
derstand. It is of course possible to integrate up and thus 
obtain the double and single differential rates, and the 
total probability. These will be the subject of a future 
study. 
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